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Abstract 



We present exact calculations of the zero-temperature partition function 
(chromatic polynomial) and W(q), the exponent of the ground-state entropy, 
for the (/-state Potts antiferromagnet with next-nearest-neighbor spin-spin 
couplings on square lattice strips, of width L y = 3 and L y = 4 vertices and 
arbitrarily great length L x vertices, with both free and periodic boundary 
conditions. The resultant values of W for a range of physical q values are 
compared with each other and with the values for the full 2D lattice. These 
results give insight into the effect of such non-nearest neighbor couplings on 
the ground state entropy. We show that the q = 2 (Ising) and q = 4 Potts 
antiferromagnets have zero-temperature critical points on the L x — > oo limits 
of the strips that we study. With the generalization of q from Z + to C, we 
determine the analytic structure of W(q) in the q plane for the various cases. 
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I. INTRODUCTION 



The g-state Potts antiferromagnet (AF) with the usual nearest-neighbor spin-spin cou- 
plings |],|[ exhibits nonzero ground state entropy, S > (without frustration) for sufficiently 
large q on a given lattice A. This is equivalent to a ground state degeneracy per site W > 1, 
since Sq = ks In W. Such nonzero ground state entropy is important as an exception to 
the third law of thermodynamics ||,|J. A physical example of nonzero ground state en- 
tropy is ice f5[- @. In this g-state Potts antiferromagnet at T = 0, the value of each spin 
must be different than the values of all of the other spins to which it is coupled. There is 
a close connection with graph theory here, since the zero-temperature partition function of 
the above-mentioned g-state Potts antiferromagnet on a lattice A or, more generally, a graph 
G, satisfies 

Z(G,q,T = Q) PAF = P(G,q) (1.1) 

where P(G, q) is the chromatic polynomial expressing the number of ways of coloring the 
vertices of the graph G with q colors such that no two adjacent vertices have the same color 
(for reviews, see |J- |T|). The minimum number of colors for which this coloring is possible, 



i.e. the minimum integer value of q for which P(G, q) is nonzero, is denoted the chromatic 
number of G, x(G). 

From eq. Ql.ll) , it follows thatQ 

W{{G},q) = lim P{G,q) 1/n (1.2) 

n— >oo 

where n = v(G) is the number of vertices of G and {G} = lim^oc G. 

Since P(G, q) is a polynomial, one can generalize q from Z + to C. The zeros of P(G, q) in 
the complex q plane are called chromatic zeros; a subset of these may form an accumulation 
set in the n — > oo limit, denoted B, which is the continuous locus of points where VK({G}, q) 
is nonanalytic. The maximal region in the complex q plane to which one can analytically 
continue the function W({G},q) from physical values where there is nonzero ground state 
entropy is denoted R\. The maximal value of q where B intersects the (positive) real axis 
is labelled q c ({G}). This value is important since W({G},q) is a real analytic solution for 



At certain special points q s (typically q s = 0, 1, .., x(G)), one has the noncommutativity of limits 
lim 9 _ >gs lim n _ >00 P(G, q) 1 ^ 71 ^ lim n _ ) . 00 \im q ^ qs P(G, q) x ^ n , and hence it is necessary to specify the order 
of the limits in the definition of W({G}, q s ) |l2|. We use the first order of limits here; this has the advantage 
of removing certain isolated discontinuities in W . 



1 



real q down to q c ({G}). For regions other than Ri, one can only determine the magnitude 



^({Ct}, q)\ unambiguously |12|. In addition to ||- [jl2[| , some previous works on chromatic 
polynomials include [[y|- [j45fl . 

In previous works we have carried out comparative studies of W for different lattices and 



have explored the effects of different lattice properties such as coordination numbers [12 



2"6fl - fl26| , |29|] - |3i| , p^ , |35| , |39[| - f4"2"fl , |45| . In general it was found that as one increased the 



lattice coordination number, the ground state entropy of the g-state Potts antiferromagnet (if 
nonzero for the given value of q), decreased. This can be understood as a consequence of the 
fact that as one increases the lattice coordination number, one is increasing the constraints 
on the coloring of a given vertex subject to the constraint that other vertices of the lattice 
adjacent to this one (i.e. connected with a bond or edge of the lattice) have different colors. 
Another way in which to explore this effect is to consider non-nearest-neighbor spin-spin 
couplings. Again, in general, these increase the constraints on the values that any given 
spin can take on, and hence decrease the ground state entropy. We wish to make this more 
quantitative and shall do so here using exact solutions for the chromatic polynomials and 
rigorous bounds. A natural starting point for studies of such non-nearest-neighbor spin-spin 
couplings is to consider the model on a given lattice and add next-nearest-neighbor (nnn) 
spin-spin couplings. Equivalently, we can redefine the lattice itself by considering it as a 
graph G with vertices at the usual lattice sites but with bonds (= edges in graph theory 
nomenclature) consisting not only of the usual bonds joining these lattice sites but also 
bonds joining next-nearest-neighbor lattice sites. We then consider the nearest-neighbor 
Potts antiferromagnet on this redefined lattice. 

Perhaps the simplest case that one could consider is the g-state Potts antiferromagnet 
in one dimension; the lattice here is just the a line T n or circle C n for the case of free and 
periodic boundary conditions, respectively (denoted FBC X and PBC X ). One has P(T n , q) = 
q(q — l) n ~ l so W = q — I and Ri is the entire q plane. For the circle, P(C n , q) = (q — l) n + 
(q — 1)(— I ) n . In this case, if \q— 1| > 1 then W = q — 1, while if \q — 1| < 1, then \W\ = 1, 
so that q c = 2. Hence for either FBC X or PBC X there is nonzero ground state entropy 
So = fcfiln^ — 1) for q > 2. The addition of next-nearest neighbor bonds converts T n or 
C n to a open or cyclic strip of triangles, respectively, with each pair sharing an edge. We 
denote these strips tri(L y = 2), N t ,BC x where BC X = FBC X or PBC X . In the cyclic case, 
the degree A (number of neighboring vertices) of each vertex is changed from 2 to 4, and 
this is also true of the internal vertices in the open case. The chromatic numbers in these 
cases are (i) % = 2 for the line T n ; (ii) % = 2 (3) for C n with n even (odd); (iii) x — 3 for 
the open triangular strip; and (iv) \ — 3 (4) for the cyclic triangular strip with the number 
of triangles N t even (odd). For the Potts antiferromagnet with nnn spin-spin couplings on 
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the line ( equivalent ly, on the open triangular strip) with n vertices, 

P(sq d (L y = 1), FBC y , FBC X , q) = P(trt(L y = 2), FBC y , FBC X , q) = q(q - l)(q - 2) n ~ 2 

(1.3) 

and hence W = q — 2, and Ri is the full q plane. 

For the Potts antiferromagnet with nnn spin-spin couplings on the circuit, the equivalence 
is with the cyclic triangular strip: 



P{sq d {L y = 1), FBCy, PBC X , q) = P{tn{L y = 2), FBC y , PBC X , q) 



;i.4) 



If the strip length involves an even number of triangles N t — n — 2m, then [35,38 



P(tri(L y = 2),N t = 2m, FBC y , PBC X , q) = q 2 - 3q + 1 + (q - 2) 2m + (q - 1) (A t2 , 3 ) m + (A t2j4 

(1.5) 



while for odd N t = n = 2m + 1 [IS] 

P(trt(L y = 2),N t = 2m + l, FBC y , PBC X , q) = -(q 2 - 3q + 1) + (q - 2) [(q - 2) 2 } m - 



_(g_l)( g _ 3 ) 



where 



(A t2 , 3 ) m + (A* 2 , 4 ) m + 



(At2,3) m — (At2,4) T 



A+2.3 — A 



t2A 



A 



42,(3,4) 



5 - 2q ± v/9-4g 



;i.6) 



;i-7) 



In both cases, g c = 3 and W = q — 2 for g > 3 (and more generally, for q in the region Ri 
given in p5|). Thus, the addition of next-nearest-neighbor couplings increases the value of 
g beyond which there is nonzero ground state entropy from 2 to 3 and decreases the value 
of the resultant entropy from So = ks ln(g — 1) to So = ks ln(g — 2) for q > 3. 

We proceed to consider the Potts antiferromagnet on the square lattice, and again add 
next-nearest-neighbor couplings, or equivalently redefine the lattice so that the bonds con- 
sists not just of the usual horizontal and vertical bonds, but also of bonds connecting the 



diagonally opposite vertices of each square. Following our earlier notation we shall 
denote this lattice as sqj, where the d refers to the addition of these diagonal bonds. For 
the square (sq) and sqd lattices, the chromatic numbers are 



X(sq) = 2 , x(sq d ) = 4 . 



'U 
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No exact solution is known for W(q) on the sq d lattice. In the absence of such an exact 
solution, Tsai and one of us (RS) have carried out Monte Carlo measurements of W(sq d , q) 
and have derived a rigorous lower bound [|2"5| 



W(sq d ,q) > 



(g-2)(g-3) 
q-1 



for q > x(sq d ) ■ 



1.9) 



This lower bound was compared with the actual value of W(sq d ,q), as determined by the 



Monte Carlo measurements for 5 < q < 10 [0,0 and was found to lie very close it (cf. 
Table III of [23]). For example, for q = 6 and q — 8, the ratio of the lower bound divided 
by the actual value was 0.981 and 0.995, and it increased monotonically toward unity as q 
increased. Since it is possible to obtain exact analytic solutions for W on infinite-length, 

one has an alternate 



finite-width strips of 2D lattices @, |2|, ||- g|, 
way to investigate W(sq d ,q), namely to calculate W exactly on strips of the sq d lattice, 
with various boundary conditions. It has, indeed, been found |J]J that for the square and 
triangular lattices, the values of W for such infinite-length strips of even rather modest 
widths are close to the corresponding values for the 2D thermodynamic limit, for moderate 
values of q. In the present work we report exact calculations of P(q), W(q), and B on strips of 
the sq d lattice with various boundary conditions. The longitudinal and transverse directions 
on the strip are taken to be x and y, respectively. In Fig. 1 we show some illustrative strips 
of the sq d lattice. 
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FIGURES 



FIG. 1. Illustrative strip graphs of the sqd lattice: (a,b), L y = 2,3 (FBC y , PBC X ) (cyclic); (c) L y = 3, 
(PBCy, PBC X ) (toroidal boundary conditions). 

An important property is that with the two added diagonal bonds, each square of the 
sqd lattice constitutes a complete graph on four vertices. (Here, the complete graph on r 
vertices, K r , is defined as the graph each of whose vertices is connected to all of the other 
r — 1 vertices by bonds (= edges); it has chromatic number x(Kr) = r -) Compared with the 
square lattice, for a given g-coloring of the sqd lattice, the addition of these bonds clearly 
increases the constraints on the coloring of each vertex and therefore decreases P(G, q). As 



we proved earlier [26|, if a lattice A' can be obtained from another, A, by connecting disjoint 
vertices of A with bonds, then W(A', q) < W(A, q) for g-colorings of the two lattices. An 
example of the application of this theorem was given in [p6fl : for q colorings of the square, 
triangular (tri), and honeycomb (he) lattices, W(tri, q) < W(sq, q) < W(hc,q). (W(sq,q) is 
strictly less than W(hc, q) except at the value q — 2, where W(sq,2) = W(hc,2) = 1.) In 
the present context, we note the inequality for g-colorings of these lattices: 

W(sq d , q) < W(tri, q) < W{sq, q) < W{hc, q) (1.10) 

(for q values where such colorings are possible). 

We use the symbols FBC^ and PBC y for free and periodic transverse boundary conditions 
and, as above, FBC X , PBC X , and TPBC X for free, periodic, and twisted periodic longitudinal 
boundary conditions. The term "twisted" means that the longitudinal ends of the strip 
are identified with reversed orientation. These strip graphs can be embedded on surfaces 
with the following topologies: (i) (FBC^FBC^): open strip; (ii) (PBC^FBCx.): cylindrical; 
(iii) (FBC y ,PBC x ): cylindrical (denoted cyclic here); (iv) (FBCj,,TPBC x ): Mdbius; (v) 
(PBC y ,PBC x ): torus; and (vi) (PBCj,,TPBC x ): Klein bottle.0 

The labelling of the strips generally follows our earlier labelling conventions. Thus for a 
strip with free transverse and longitudinal boundary conditions, (FBC y , FBC X ), the length 
of the strip is taken to be m+1 squares or equivalently edges, with L x = m+2, and the width 
is L y vertices. This strip thus has n = L x L y vertices and e = 4:L x L y — 3(L X + L y ) + 2 edges. 
For cyclic strips, the width is defined in the same manner and the length is L x vertices 
or equivalently edges. For strips with periodic transverse boundary conditions, including 



These BC's can all be implemented in a manner that is uniform in the length L x ; the case (vii) 
(TPBC y ,TPBCa;) with the topology of the projective plane requires different identifications as L x varies 
and will not be considered here. 
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(PBCy, FBC X ) and (PBC y , PBC X ), a width of L y = 3 means that the cross section involves 
L y vertices. For L x = m > 3 to avoid certain degenerate cases, the sq^ strips with either 
cyclic or torus boundary conditions have n = L x L y vertices. With the same restriction, the 
cyclic strips have e = L x (AL y — 3) edges and the torus strips have e = 4L x L y edges. 

Let us next comment on the planarity or nonplanarity of the strips of the sqd lattice with 
various boundary conditions. We shall concentrate here on nondegenerate cases where the 
strips are proper graphs without multiple edges. Consider first the strips with (FBC y , FBC X ) 
boundary conditions. For L y = 2 and arbitrary L x , it is easy to show that these are planar 
by taking the second diagonal bond for each square and drawing it external to the strip; by 
redefining the labelling of the x and y axes, it follows that this strip with (FBC y , FBC X ) 
boundary conditions, L x = 2, and arbitrary L y is also planar. For other cases we shall make 
use of two theorems from graph theory. The first of these states that if G is a planar graph 
with n vertices and e edges with n > 3, then e < 3(n — 2) (e.g. Corollary 11.1(c) in |lQfl) and 



the second states that if G is a planar graph with n > 4, then G has at least four vertices of 
degree A < 5 (e.g. Corollary 11.1(e) in |10|). Now 



3(n - 2) - e = -L x L y + 3{L X + L y ) - 8 for sq d , (FBC y ,FBC x ) (1.11) 

so that for sufficiently great L x and/or L y , 3(n — 2) — e is negative and hence the strip is 
nonplanar, by the first theorem. For example, 3(n — 2) — e < if L y = 4 and L x > 5 or vice 
versa, i.e., L x = 4 and L y > 5; and similarly, if L x = L y > 5. Considering next the strips of 
the sqd lattice that are cyclic, i.e., have (FBC y , PBC X ) boundary conditions, we have 

3 ( n _ 2 ) - e = -L x L y + 3L x -6 for sq d , (FBC y , PBC X ) . (1.12) 

Hence, 3(n — 2) — e < for all L x if L y > 3, so that these strips are nonplanar. For the 
strips with (PBC y , PBC X ), i.e., torus boundary conditions, we have 

3( n -2)-e = -(L x L y + 6) for sq d , (PBC y ,PBC x ) (1.13) 

so that these strips are also nonplanar. This can be seen alternatively by observing that each 
vertex on the torus strips has degree A = 8 and applying the second theorem cited above. 
The second theorem also shows that the sqd strip with Klein bottle boundary conditions is 
nonplanar. 

A generic form for chromatic polynomials for recursively defined families of graphs, of 
which strip graphs G s are special cases, is 

p((G s ) m , q ) = J2 cg s ,M(*g s MT (i-w) 

3=1 
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where ca at j(q) and the Ng„,x terms \G s ,j(q) depend on the type of strip graph G s , as indicated, 
but are independent of m. 



II. STRIPS WITH (FBC y ,FBC XJ 



A. L y = 2 



The chromatic polynomial for the strip of the sqd lattice with L y = 1 and free transverse 
and longitudinal boundary conditions was given above in (|1.3|) . For the L y = 2 case the 
chromatic polynomial is 



P(sq d (L y = 2) m , FBCy, FBC X , q) = q{q - 1) \{q -2){q- 3) 
In the m — > oo limit, 



m+l 



W(sq d (L y = 2), FBCy, FBC X , q) = {q - 2){q - 3) 



1/2 



(2.1) 



(2.2) 



with B 



B. Ly — 3 



For the L y = 3 strip, we use the same generating function method as we have before 
1 29| , |30|J35| . In general, for the family of strip graphs G s , the generating function T(G s ,q,x) 



is a rational function of the form 



with 



r(G„^) = ^4 (2.3) 



Af(G s ,q,x) = J2 A G*>^ xj ( 2 - 4 ) 

j=0 



and 



V(G s ,q,x) = l + J2b Ga ,Mx J (2-5) 

3=1 

where the Ac Bt i and bc s ,« are polynomials in q (with no common factors) and 

d N = deg x {M) (2.6) 
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and 



d v = deg x (V) . 



(2.7) 



This generating function yields the chromatic polynomials via the Taylor series expansion in 
the auxiliary variable x: 



T(G s ,q,x) = J2 P (( G s)m,q)x r 



{21 



m=0 



where we follow the notational convention in |[29|| , according to which a strip is considered to 
be comprised of m repetitions of a basic subgraph unit H connected to an initial subgraph 
/; here we take I = H so that a strip with a given value of m has m + 1 columns of K^s and 
m + 2 vertices in the longitudinal direction. The denominator can be written in factorized 
form as 



V Gs = - X Gsd x) . 



(2.9) 



These are the Ag s j's in eq. ( |1.14|) ; the coefficients are determined by eqs. (2.14) or (2.19) in 
For L y = 3, we find dx> = 2, dj^ = 1, and 



T{sq d {L y = 3), FBC y , FBC X , q, x) 



q(q -l)(q- 2)(q - 3) 2 (q - 2) - (q - l)(q - 3)x 



1 - 



3)(g 2 - 6q + U)x + (q-2)(q- 3) 3 x 2 



(2.10) 



The denominator can be written as 

T^sqdZo — (1 — ^sqd3o,lX)(l ~ \ S qd3o,2%) 

where 

A S9d30)(1) 2) = \{q ~ 3) [q 2 - Qq + 11 ± (q 4 - 12g 3 + 54g 2 - 112g + 97^ 

and the shorthand sqd3o denotes the strip of the sqd lattice with L y = 3 and open (o) bound- 
ary conditions. From this, using the general formulas in [3C], one can write the chromatic 
polynomial in the form of eq. ( |1.14|) (with N\ = 2). In the m — > oo limit, 



1/2-, 



(2.11) 



(2.12) 



W{sq d {L y = 3),FBC y ,FBC x ,q) = (X sqd3o> 



,1/3 



(2.13) 
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FIG. 2. Locus B for W for the 3 x oo strip of the sqd lattice with free transverse and longitudinal 
boundary conditions. Chromatic zeros are shown for the case L x = 20 (i.e., n = 60). 

The nonanalytic locus B is shown in Fig. |^ and is comprised of an arc stretching between 
endpoints at q ~ 1.95 + 1.432 and 4.05 + 0.396i, together with the complex conjugate arc. 
In agreement with the general discussion given before [29,30,¥l] , these four points are the 
branch points of the square root in eq. ( p,12| ). B does not intersect the real q axis, so that 
no q c is defined. The region Ri is the entire q plane, with the exception of the arcs lying on 
B. 
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Here we find dx> = 4, d_\f = 3. Again using the shorthand notation sqdAo to denote this 
open strip of the sqd lattice with L y = 4 we have 

b sqd 4o,i = -<? 4 + 13g 3 - 68q 2 + lUq - 176 (2.14) 
b sqd ± ,2 = (q- 3)(2g 5 - 33g 4 + 219g 3 - 729g 2 + 1214g - 803) (2.15) 
b sg Mo,3 = (q~ 3) 3 (g 5 - 17g 4 + 118g 3 - 420g 2 + 770g - 586) (2.16) 

bsgdAoA = ~(q ~2){q- 3) 6 (g - 4) . (2.17) 

For the functions A sq d4, ,j in the numerator, Af, it is convenient to extract a common factor 
and thus define 



Then 



A sqdA o,j = q{q -l)(q- 2)(g - 3) 3 A sqMoJ . (2.18) 

A sqdAOt0 = {q- 2) 2 (2.19) 

A sqMO)l = -(2g 4 - 21g 3 + 78g 2 - 113g + 47) (2.20) 

A sqMo , 2 = ~{q ~ 3)(g 5 - 14g 4 + 77g 3 - 204g 2 + 245g - 91) (2.21) 



and 



A sqdAo , 3 = {q- l) 2 {q - 3) 3 (g - 4) . (2.22) 
Let us write the denominator as 

4 

T^sqdAo = ~ KqdiojX) ■ (2.23) 

Then 

W = (A S3d40)i , m ax) 1/4 for qeR-y (2.24) 



where \ S q d 4o,j,max is the X sq d4o,j m (|2.23|) with maximal magnitude in this region. 
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FIG. 3. Chromatic zeros for the 4 x L x strip of the sqd lattice with free boundary conditions and L x = 20 
(i.e., n = 80). 

Chromatic zeros are shown in Fig. [5] for L x = 20, i.e., n = 80. For this great a length, 
these chromatic zeros give a reasonably good approximation to the asymptotic locus B. As 
is evident from this figure, B does not cross the real q axis, so that no q c is defined. 

III. STRIPS WITH (FBCy, (T)PBC X ) 

A. L y = 2 

The chromatic polynomial for the L y — 1 cyclic strip of the sqd lattice was given above in 
eqs. ( |1.5| ) and Here we consider the L y = 2 strip of the sqd lattice with (FBC y , PBC X ), 
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i.e. cyclic, boundary conditions. For a given value of L x , this cyclic strip graph is identical 
to the corresponding strip with Mobius boundary conditions (FBC y ,TPBC x ): 

G(sq d , L y = 2, L x , FBCy, PBC X ) = G(sq d , L y = 2, L x , FBC y , TPBC X ) (3.1) 

This can be proved by calculating the adjacency matrices for the cyclic and Mobius strips, 
which are identical. (Here the adjacency matrix of an n-vertex graph is the n x n ma- 
trix A with Ay equal to the number of bonds (edges) that connect the i'th and j'th 
vertices. The adjacency matrix fully defines the graph.) Because of the identity of 
the L y = 2 cyclic and Mobius strips, we shall refer to them both with the designation 
sqd(L y = 2). m , FBCy, (T)PBC X ). For the general cyclic strip of the sqd strip, with L x > 4 
to avoid degenerate cases, the chromatic number is given by 

X {^L t ,L z ,FBC m PBC x ) = [\ H'X^i (3 ' 2) 

For the present L y = 2 cyclic/Mobius strips, the degenerate cases are as follows: for L x = 2, 
the strip reduces to K A while for L x — 3 it reduces to K 6 , with x(Kp) = P- 

The chromatic polynomial for the cyclic strip of the sqd lattice with L y = 2 and L x = in 

is 



20 



P(sq d (L y = 2), FBCy, (T)PBC x ) m , q) = ~q(q - 3)2 m + [(q - 2)(q - 3)] m + (q - 1)[2(3 - q)\ 

(3.3) 

We determine the boundary B to be the union of a circle centered at q = 2 with radius 2 
and a circle centered at q = 3 with radius 1: 

B:{\q-2\ = 2} U {|g-3| = l}. (3.4) 

These two circles osculate (i.e., intersect with equal tangents) at q c , where 

q c (sq d (Ly = 2), FBC y , (T)PBC X ) = 4 . (3.5) 

In the terminology of algebraic geometry, this point q c is thus a tacnode. 
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FIG. 4. Locus B for W for the 2 x oo cyclic or Mobius strip of the sqa lattice. Chromatic zeros are 
shown for the case L x = 20 (i.e., n = 40). 

This locus is shown in Fig. |j. The locus B separates the q plane into three regions: (i) 
the outermost region Ri which is the exterior of the larger circle, \q — 2| > 2 and which thus 
includes the real intervals q > 4 and q < 0; (ii) region R2 which is the interior of the smaller 
circle, I*? — 3| < 1 and includes the real interval 2 < q < 4; and (iii) region R3 which is the 
interior of the larger circle | q — 2 1 < 2 minus the smaller disk | q — 3 1 = 1 and includes the real 
interval < q < 2. Thus, B crosses the real q axis at q = 0, 2, 4 and q c = 4. As is evident in 
Fig. |], the chromatic zeros lie near to the asymptotic locus B. In the various regions 

W = [(g-2)(g-3)] 1/2 for q e Ri (3.6) 
\W\ = 2 1/2 for q G R 2 (3.7) 
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\W\ = |2(g-3)| 1/2 for q G R 3 (3.8) 

(for q in regions other than Ri, only the magnitude can be determined unambigu- 
ously). 

We define the sum of the coefficients as 

C(G) = J2cg,j. (3.9) 



For sufficiently large positive integer q, the coefficient cqj in ( |3.9p can be interpreted as the 
multiplicity of the corresponding eigenvalue Xgj of the coloring matrix, i.e., the dimension 
of the corresponding invariant subspace in the full space of coloring configurations |IR 



We recall that the coloring matrix can be defined as the matrix whose i,j element is 1 (0) 
if the coloring configurations on two adjacent transverse slices of the strip are compatible 
(incompatible). Thus, in the absence of zero eigenvalues of the coloring matrix, C(G) is the 
dimension of the space of coloring configurations of such a transverse slice. For the cyclic 
strip graphs of the sqd lattice, the transverse slice is the line graph L n of length L y vertices, 
so the space of coloring configurations of this transverse slice is 

P(L n ,q) = P(T n ,q) = q(q-ir . (3.10) 

In cases where the cyclic strip and the Mobius strip of the sqd lattice are identical, the full 
coloring matrix automatically takes account of both contributions, so that for each individual 
strip, corresponding to a given permutation in the identification of vertices at the longitudinal 
boundary, one must divide by the symmetry factor. In the present case, L y = 2 and there 
are two permutations of the identifications of the boundary conditions that give identical 
strip graphs, so that this symmetry factor is 1/2! so that 

C(sq d , L y = 2, FBCy, PBC X ) = C(sq d , L y = 2, FBC y , TPBC X ) = l -q(q - 1) (3.11) 

This agrees with the sum of the coefficients in the expression (|3.3|) . A remark that will be 
relevant later is that if the coloring matrix has a zero eigenvalue of multiplicity c zero , then, 
since this eigenvalue does not appear in (|1.14j) , the sum of the coefficients that do appear 
in the chromatic polynomial (|1.14|) is equal to the full dimension of the space of coloring 
configurations minus c zero . 
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B. Ly — 3 



We have calculated the chromatic polynomial for the next wider cyclic strip, with L y = 3. 
For this strip the chromatic numbers \ are the same as for the L y = 2 sq d strip. We find 

N sqd ,L y =3,cyc.,x = 16 and 
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P(sq d (L y = 3), FBCy, PBC x ) m , q) = ^ c sqd ^(X sqd ^) m (3.12) 

3=1 

(The terms X S qd3c,j are the same for cyclic and Mobius longitudinal boundary conditions.) 
With the \ sqd 3c,/s ordered according to decreasing degrees of their coefficients for the cyclic 
strip, we find 

A S gd3c,l = 1 (3.13) 

K q d3c,2 = —2 (3-14) 

Kqd3c,3 = — 3 (3.15) 

A S gd 3c ,4 = g-3 (3.16) 



Kqd3c,5 = -(q - 3) (3.17) 



and 



Kqd3c,( 6 ,7) = q - 4 ± ^2qi - Uq + 25 . (3.18) 
The terms \ sqd 3c,j, j = 8, 9, 10 are the roots of the cubic equation 

e - 4(g - A)^ 2 + (q 2 - lOq + 17)^ + 2(g - l)(g - 3) . (3.19) 
For j = 11 we have 

A^scii = -(? " 3) 2 • (3.20) 
The terms \ sqd 3 C ,j, j = 12, 13, 14 are the roots of the cubic equation 

e + (2q 2 - 15q + 30)£ 2 - (q - 3) 2 (g 2 - 5q + 5)£ - 2(g - 3) 4 . (3.21) 



16 



Finally, 

Kqd3c,(15,16) = ^sqd3o,(l,2) ■ (3.22) 

For the coefficients we calculate 

c agd 3ci = i(g-l)(g-2)(g-3) (3.23) 

c sgd 3 C>2 =^q(q-2)(q-A) (3.24) 

c sqd 3c,3 = \q{q ~ l)(g - 5) (3.25) 
o 

c^ 3c ,i = ^g(g - 3) f o r J =4, 8, 9, 10 (3.26) 

c sg d3c,j = l)(?-2) for J = 5, 6, 7 (3.27) 



c s9d3Cli = g-l for j = 11, 12, 13, 14 (3.28) 



and 



c sq d3c,j = 1 for j = 15, 16 . (3.29) 
Summing the coefficients c sq d3 C ,j, we find 

C(sq d (L y = 3), FBCy, PBC X ) = ^q(q + l)(4q - 7) . (3.30) 

Since for L y = 3 the cyclic and Mobius strips of the sqd lattice are distinct, the full sum 
of eigenvalue multiplicities is equal to ( |3.10|) with n = L y = 3 (without dividing by any 



symmetry factor). The sum of the coefficients appearing in ( |3.12j ) is less than this quantity, 
q(q — l) 2 , by the amount 

C S qd3c,zero = ^?(2? 2 - 9q + 13) (3.31) 

6 

which indicates that for this strip the coloring matrix has a zero eigenvalue with this multi- 
plicity. 
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FIG. 5. Locus B for W for the 3 x oo cyclic or Mobius strip of the sqa lattice. Chromatic zeros are 
shown for the case L x = 20 (i.e., n = 60). 

The boundary B is shown in Fig. [5]. This boundary separates the q plane into four 
regions. These include (i) the outermost region R±, which contains the semi-infinite intervals 
q > q c and q < 0, where q c is 

q c {sq d {L y = 3), FBC y , PBC X ) = 4.254654... (3.32) 

(a root of the equation q 4 — 14g 3 + 72q 2 — 168g + 162 = 0); (ii) a narrow crescent-shaped 
region R2 containing the real interval 4 < q < q c \ (iii) the region R3 containing the real 
interval 2 < q < 4; and (iv) the region R^ containing the real interval < q < 2. Associated 
with these regions are two complex-conjugate pairs of triple points, as is evident in Fig. [5[ 
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Note that q c is not a tacnode for the (L x — > oo limit of the) L y = 3 cyclic strip, in contrast 
to the situation for the corresponding L y = 2 cyclic strip. 
In the various regions 

W = (\ sqd3c , 15 ) 1/3 for qeR! (3.33) 
\W\ = 3 1/3 for q e R 2 (3.34) 

and 

\W\ = |A sgd3c ,8io m r /3 for q e R 3 (3.35) 
where \ S qd3c,8Wm denotes the root of the cubic ( |3.19 ) of maximal magnitude in R 3 , and 



\W\ = \X sqd3c ,ni4m\ 1/3 for q G R A (3.36) 
where A S(?( i3 Ci i2i4m denotes the root of the cubic ( |3.21|) of maximal magnitude in i? 4 . 



IV. STRIPS WITH (PBC y ,FBC x ) 
A. L y = 3 

For the L y = 3, L x = m + 2 strip of K^s forming squares, with (PBC^FBGr) boundary 
conditions (for which n = 3(m + 2)) we find 

-| m+l 

P(sq d (L y = 3) m , PBCy, FBC X , q) = q(q -!)(<?- 2) (q - 3)(g - 4)(g - 5) (4.1) 



whence 



W(sq d (Ly = 3),PBC yi FBC x ,q)= (q - 3)(q - 4)(g - 5) ^ . (4.2) 



The continuous nonanalytic locus B = W has isolated branch point singularities where it 
vanishes at g = 3,4, and 5. Aside from these, Ri is the full q plane. 



B. Ly = 4 



For this case it is again convenient to give our results in terms of a generating function 
T(sq d (Ly = 4), PBCy, FBC x ,q,x). We find d-p = 3, = 2, and (using the abbreviation 
sqdicyl for this strip) 
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b sqdAcy i,x = -g 4 + 16g 3 - 104g 2 + 316g - 372 (4.3) 

b sgMcy i,2 = (q~ 3)(5g 4 - 74g 3 + 422g 2 - HOOg + 1109) (4.4) 

and 

b sgdW = ~{q - 2)(g - 3) 2 (2g 2 - 16g + 33) . (4.5) 

Denning 

A aq dAcyl,j = <?(<? ~ ~ 2 )(<? - 3) A sgd 4 c2/ / ; j (4.6) 

we calculate 

A s?d4c3/ / )0 = q A - 14g 3 + 79g 2 - 210g + 220 (4.7) 

A sqMc yi,i = -(5g 5 - 79g 4 + 501g 3 - 1586g 2 + 2485g - 1513) (4.8) 

and 

A sqMc yl,2 = (q - 3) {q 2 - 3q + 3) (2g 2 - 16q + 33) . (4.9) 

Writing 

3 

T^sqdAcyl = ~~ Kqd4,cyl,j%) (4-10) 

we have 

W = (XsqdAcyl^max) 1 ^ for q E R\ (4.11) 

where X S qMcyl,j,max is the \ sq dAcyi,j hi ( |4.10|) with maximal magnitude in this region. 
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FIG. 6. Chromatic zeros for the 4 x L x cylindrical strip of the sqd lattice, i.e. with (PBC y , FBC X ) 
boundary conditions with L x — 20 (i.e., n — 80 vertices). 

Chromatic zeros for the 4 x L x strip of the sqd lattice with cylindrical boundary conditions 
are shown in Fig. ^ for L x = 20, i.e. n = 80. Again, the length L x is sufficiently great that 
these give a good idea of the location of the asymptotic curve B. Since B does not cross the 
real q axis, there is no q c for this case. 



V. STRIP WITH L y = 3 AND (PBC y , (T)PBC X ) 



We consider here the sqd strip with L y = 3 and torus boundary conditions 
(PBCy, PBC X ). By the same method as mentioned above, e.g., calculating the associated 
adjacency matrices and showing that they are the same, it follows that for a given L x , this 
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strip with torus boundary conditions is identical to the corresponding L y = 3 strip with 
Klein bottle boundary conditions (PBC y ,TPBC x ): 

G(sq d , L y = 3, L x , PBC y , PBC X ) = G(sq d , L y = 3, L x , PBC y , TPBC X ) (5.1) 

Since there are 3 vertices on the vertical slice, and hence 3! permutations that yield identical 
graphs, it follows, as explained above, that the sum of the eigenvalue multiplicities for each 
individual strip contributes only 1/3! of this total: 

C(sq d , L y = 3, PBCy, (T)PBC X ) = ip(C 3 , q) = \q{q - l)(g - 2) (5.2) 
These strips have the chromatic number 

x(sq d (L y = 3) m , PBCy, (T)PBC X ) = { * ^ (5-3) 

For m = 2 and m — 3, the L y = 3 sq d torus strip degenerates to K 6 and iTg, respectively, 
with x{Kp) = P as before; for m = 5 this strip has x — 8- The fact that the values of \ f° r 
the L y = 3 torus strip of the sq d lattice in eq. ( |5.3|) and the special cases noted are larger 
than the value \ = 4 for the infinite 2D sq d lattice can be ascribed in part to the constraints 
arising from the small girth of the triangular transverse cross section of these strips. 

We calculate the chromatic polynomial by iterated use of the deletion-contraction theo- 
rem, via a generating function approach p9| , p0| . From this we obtain the chromatic polyno- 
mial in the form (|1.14 ) via the general formulas in |30] and obtain 



P(sq d (L y = 3) m , PBCy, (T)PBC X , q) = ±q(q - l)(q - 5)(-6) m + \q{q - 3)[6(g - 5)] m 

r i m r i m 

+(q-l) _3(g-4)(g-5) + (q - 3)(q - A){q - 5) . (5.4) 



Thus, N\ = 4 for this strip. The labelling of the coefficients Cj and terms Xj in eq. (5~3) 
is consecutive. Explicitly calculating the sum of the coefficients in the chromatic polyno- 
mial ( |5.4j ), one sees that the result agrees with ( |5.2j ). It is interesting that the coefficients 
that occur in ( |5.4| ) are the same as a subset of the coefficients that occur in the chromatic 
polynomial for the L y = 3 cyclic strip. We find that the degrees in x of the numerator and 
denominator of the generating function for this strip graphs are 2 and 4, so that all of the 
A,'s in dn = 4 contribute to P. 
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FIG. 7. Locus B for W for the 3 x oo strip of the sq d lattice with (PBC Vl {T)PBC X ) = torus or Klein 
bottle boundary conditions. Chromatic zeros are shown for the case L x = 20 (i.e., n = 60). 

The nonanalytic locus (boundary) B is shown in Fig. [7] and consists of three circles that 
osculate at q c , where 

q c (sq d (L y = 3), PBC y , (T)PBC X ) = 6 (5.5) 

namely, 

B:{\q-3\ = 3} U {|g-4| = 2} U {\q - 5| = 1} . (5.6) 

Thus, as was true for the L y = 2 cyclic strip, q c is a tacnode. Evidently, B crosses the real 
axis at q — 0, 2, and 4 as well as at q c . This locus £> divides the q plane into four regions: 
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(i) the outermost region Ri, including the real intervals q > 6 and q < 0; (ii) region R2, 
the interior of the smallest circle, \q = 5| = 1, containing the real interval 4 < q < 6; (iii) 
region R3, the interior of the circle \q — 4| = 2 minus the disk \q — 5| = 1 comprising R2 and 
including the real interval 2 < g < 4; and (iv) region R4, the interior of the largest circle, 
\q — 3 1 = 3 minus the disk |<7 — 4| =2 and including the real interval < q < 2. We have 



W 



(?-3)(g-4)(g-5) 



1/3 



for q E R\ 



(5.7) 



The fact that this coincides with the 14 7 calculated for the corresponding strip with 
(PBCy, FBC X ) (see eq. (|4.2j)) is a general result p5| , |39| , ^T| . For the other regions we have 



\W\ = 6 1/3 for q e R 2 



(5.8) 



\W\ = |6(g-5)| 1/3 for q G R 3 (5.9) 

and 

|W| = |3(g-4)(g-5)| 1/3 for q E R A . (5.10) 

Evidently, for all of these strips, the locus B has support for Re(q) > 0. 

It is of interest to comment further on the q c values for (the L x — > 00 limits of) these strips 
of the sqd lattice. In our previous exact calculations of chromatic polynomials for various 
strip graphs of regular lattices and the resultant W functions for their L x — > 00 limits, it 
was found that if one uses free transverse boundary conditions and periodic longitudinal 
boundary conditions the value of q c for a given family is a non- decreasing function of L y . 
Our results for the strips of the sqd lattice exhibit the same behavior. Hence, our finding 
that q c ~ 4.25 for the (L x — > 00 limit of the) L y = 3 cyclic/Mobius strip graph suggests 
that q c for the infinite 2D sqd lattice, which we denote q c (sqd), is greater than 4.25. Note 
that we cannot use our finding that q c = 6 for the L x —* 00 limit of the L y = 3 torus/Klein 
bottle graph of the sqd lattice to suggest that q c (sqd) might be 6 because we have previously 
obtained exact solutions for W that show that q c is not, in general, a non-decreasing function 
of L y on strip graphs with periodic transverse boundary conditions [p9 , |4~5l . For example, 



from exact results, we have found that for the L x — > 00 limit of strips of the triangular 
lattice with cylindrical boundary condition (PBC y , FBC X ), q c = 4 for L y = 4 [29|, while 



q c ~ 3.28 for L y = 4 and q c ~ 3.25 for L y = 5 P5| . For the square, triangular, and sqd 
lattices, constructed, say, as the L x , L y — ► 00 limits of open rectangular sections, one has 
the chromatic numbers \ = 2, 3, and 4, respectively. Now the Potts antiferromagnet has a 



zero-temperature critical point at q = 3 on the square lattice 0, [17fl, [18[ and at q = 4 on the 
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triangular lattice Jjl| respectively (which should be independent of boundary conditions used 
in taking the thermodynamic limit), corresponding to the values q c {sq) = 3 and q c (tri) = 4. 
These results are consistent with the possibility that q c {sqd) = 5, i.e., the possibility that the 
q = 5 Potts antiferromagnet has a T = critical point on the sqd lattice. However, there is 
not a 1-1 correspondence between chromatic number and q c \ for example, the kagome lattice 
(again constructed, say, as the limit of a finite section with free transverse and longitudinal 
boundary conditions) has x = 3 like the triangular lattice, but q c = 3, in contrast to the 
q c = 4 value for the triangular lattice. If, indeed, q c (sqd) = 5, this would also mean that q c for 
an infinite-length, finite-width strip could be larger than q c for the full infinite lattice since 
we have obtained q c = 6 in eq. (|5.5|) from our exact solution for the chromatic polynomial 
for the strip with torus boundary conditions above. 

VI. RIGOROUS LOWER BOUNDS ON W AND APPROACH TO THE 

INFINITE- WIDTH LIMIT 

Here we present rigorous lower bounds on W for strip graphs and show that these are 
very close to the actual values obtained from our exact solutions and hence serve as very 
good approximations to the actual W functions. Using our exact solutions and these ap- 
proximations, we then determine how, for a given q, the values of W for the infinite-length, 
finite-width strips approach the value for the infinite 2D sqd lattice as the strip width L y 
gets large. 

As discussed in [^,^1, a general result for a given type of strip graph G s is 
W(G s (L y ), BC y , FBC X , q) = W(G s (L y ), BC y , PBC X , q) for q > q c (G s (L y ), BC y , PBC X ) 

(6.1) 

where BC y = FBC y or PBC y . Hence, for example, W{sq d {L y = 2), FBC y , FBC X , q) = 
W(sq d (L y = 2), FBCy, PBC X , q) for q > 4 and 

W(sqd(Ly = 3),PBCy,FBC x ,q) = W(sq d (Ly = 3),PBCy,PBC x ,q) for q>6 (6.2) 

We recall that, using coloring matrix methods flT5 |, Tsai and one of us (RS) previously 



derived rigorous upper and lower bounds on W for various 2D lattices [23]- ||26|| . It was found 
that the lower bounds Wu,, were actually very good approximations to the actual W values, 
as determined, e.g., by Monte Carlo simulations. This was also seen analytically from the 
property that the large-g Taylor series expansions of q~ x W and q~ x Wt,b. coincided to several 
orders beyond the first term, which is unity. The lower bound (|1.9|) was derived as part 
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of this work. As noted above, this bound agrees very well with the actual values of W, as 
determined via Monte Carlo measurements (see Table III of ||25|| ). 

Using the same methods, we can obtain a lower bound on W for the sqd strips of interest 
here. Here we restrict to q > 6. For the case of free transverse boundary conditions and 
either free or periodic longitudinal boundary conditions we obtain the lower bound 

i i 

W(sq d (L y ), FBC yj BC XJ q) > [fa ~ 2 )(g ~ 3 )[ . (6 . 3) 

(q - 1) L y 

For the sqd strips with periodic transverse boundary conditions and either free or periodic 
longitudinal boundary conditions, we obtain the lower bound 

W(sq d (L y ), PBC y , BC X , q) > (6.4) 

where 

= (g/gXg - 3)2^ + [(g - 2)(q - 3)] L y + (q- 1)[2(3 - q)] L « 

(q - l) L v + (q - l)(-l) L v ' 1 ' ' 



Evidently, for large g, the lower bound ( |6.4| ) with ( |6.5| ) goes over to ( |1.9|) for the infinite 
lattice. 

In Table Q, we compare the values of W for 6 < q < 10 from exact solutions for the 
L y = 2, 3, 4 open strips with the respective lower bound ( |1.9| ). For each pair of values of L y 
and q, the upper entry is the value of W from the exact solution and the lower entry is the 
ratio 

W(sq d (L y ), FBC y , BC X) q) 
WF W(sq d (L y ),FBC y ,BC x ,q) lb . 1 ' J 

where W(sqd(L y ), FBC y , BC X , q)tb. is the lower bound (£.b.) given by the right-hand side of 
eq. (|6.3|) and the numerator and denominator are independent of the boundary conditions in 
the longitudinal direction. The bound is identical to the exact expression for W for L y = 2 
and is very close for L y = 3 and L y = 4. Thus, just as was found in our earlier work with 



Tsai [24|- 1 26(1 , the lower bound is not just a bound but also a very accurate approximation 
to the exact value of W, especially for moderate and large q. Having confirmed this again 
for the present type of strip graphs, we use this approximation to study the approach to the 
infinite-width limit, i.e. the full infinite 2D sqd lattice. In Table |lj we show values of W for 
6 < q < 10 and widths 2 < L y < 10, as compared with the values for L y = oo, i.e. for the 
full 2D sqd lattice. For each pair of values L y and q, the upper entry is the value of W and 
the lower entry is the ratio of this value divided by the corresponding value of W for the 
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2D sqd lattice. For L y = 2, 3, 4, the values of W are from evaluations of the exact solutions, 
while for L y from 5 to 10 they are from the approximation provided by ( |6.4|) with ( |6.5|) and 
for L y = oo they are from the Monte Carlo measurements given in ||25|| . One sees that for 
a fixed L y , the value of W on the infinite- length strip approaches that for the 2D lattice as 
q increases. Clearly, for fixed q, the value of W for the infinite-length strip approaches the 
value for the 2D lattice as L y increases, and Table |I| provides a quantitative picture of this 
approach; for example, for L y = 10 and a moderate value of q such as 7 or 8, the W values 
are within several per cent of their 2D lattice values. We recall that this approach for this 



type of strip with free transverse boundary conditions was proved to be monotonic in [31 
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TABLES 

TABLE I. Values of W for infinite- length, finite- width strips of the sqd lattice with free transverse 
boundary conditions, as functions of q, from evaluation of our exact solutions. For each pair of values of L y 
and q, the upper entry is the value of W from the exact solution and the lower entry is the ratio Rwf ■ 



L y i q -> 


6 


7 


8 


9 


10 


2 


3.464 


4.472 


5.477 


6.481 


7.483 




1 


1 


1 


1 


1 


3 


3.083 


4.066 


5.055 


6.047 


7.0415 




1.006 


1.003 


1.002 


1.001 


1.001 


4 


2.909 


3.878 


4.857 


5.842 


6.831 




1.009 


1.004 


1.002 


1.0015 


1.001 
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TABLE II. Values of W for infinite-length, width L y strips of the sqd lattice with (FBC y ,BC x ), as 
functions of q, compared with the corresponding values for the infinite 2D sqd lattice. For each value of L y 
and the upper entry is the value of W and the lower entry is the ratio of this value divided by W(sqd, q) 
for the 2D sqd lattice. For L y = 2, 3, 4, the values of W are taken from the exact solutions; for 5 < L y < 10 
the values are from (6.3). 



L y [ q — 


6 


7 


8 


9 


10 


2 


3.46 


4.47 


5.48 


6.48 


7.48 




1.42 


1.33 


1.27 


1.23 


1.20 


3 


3.08 


4.07 


5.055 


6.05 


7.04 




1.26 


1.21 


1.17 


1.15 


1.13 


4 


2.91 


3.88 


4.86 


5.84 


6.83 




1.19 


1.15 


1.13 


1.11 


1.10 


5 


2.78 


3.75 


4.73 


5.71 


6.70 




1.14 


1.11 


1.10 


1.08 


1.07 





1. 1 1 


O.OO 


4.00 


O.OO 


b.oz 




1.11 


1.09 


1.08 


1.07 


1.06 


7 


2.665 


3.625 


4.60 


5.58 


6.56 




1.09 


1.08 


1.07 


1.06 


1.05 


8 


2.63 


3.59 


4.56 


5.53 


6.52 




1.075 


1.07 


1.06 


1.05 


1.04 


9 


2.60 


3.56 


4.53 


5.50 


6.48 




1.06 


1.06 


1.05 


1.04 


1.04 


10 


2.58 


3.535 


4.50 


5.48 


6.46 




1.06 


1.05 


1.04 


1.04 


1.04 


oo 


2.45 


3.37 


4.31 


5.27 


6.24 




1 


1 


1 


1 


1 
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We next study the approach of the values of W on strips of the sqd lattice with periodic 
transverse boundary conditions to their values for the infinite 2D sqd lattice. In Table |T|, we 
compare the values of W for 6 < q < 10 from our exact solutions for the L y = 3, 4 cylindrical 
strips with the respective lower bound (|6.4j) with ( |6.5|) . For each pair of values of L y and q, 
the upper entry is the value of W from the exact solution and the lower entry is the ratio 

W(aq d (L y ),PBC yt BC x ,q) 
WP W(sq d (L y ),PBC y ,BC x ,q) t . b . 1 ' ' 

where W(sqd(L y ), PBC y , BC X , q)tb. is the lower bound (£.b.) given by the right-hand side of 
eq. ( |6.4| ) with ( |6.5| ) and the numerator and denominator are independent of the boundary 
conditions in the longitudinal direction. The bound is identical to the exact expression for 
W for L y = 3 and is very close for L y = 4. Hence, as was the case for the open strip of the 
sqd lattice, the lower bound is not just a bound but also a very accurate approximation to 
the exact value of W, especially for moderate and large q. Having established this, we use 
this approximation to study the approach to the infinite-width limit, i.e. the full infinite 2D 
sq d lattice. In Table |IV| we show values of W for 6 < q < 10 and widths 3 < L y < 10, as 
compared with the values for L y = oo, i.e. for the full 2D sqd lattice. For each pair of values 
L y and q, the upper entry is the value of W for the infinite cylindrical or torus strip and the 
lower entry is the ratio of this value divided by the corresponding value of W for the 2D sqd 
lattice. For L y = 3, 4, the values of W are from evaluations of the exact solutions, while for 
L y from 5 to 10 they are from the approximation provided by bound ( |5.3| ) and for L y = oo 
they are from the Monte Carlo measurements given in [p5| . One sees that for a fixed L y , the 
value of W on the infinite-length strip approaches that for the 2D lattice as q increases. Also, 
for fixed q, the value of W for the infinite-length strip approaches the value for the 2D lattice 
as L y increases. In |[l| it was shown that this approach is non-monotonic for strips of the 
square and triangular lattice with periodic transverse boundary conditions and the same is 
true here. The approach of the values of W for the infinite-length finite-width strips to their 
respective values for the infinite sqd lattice is more rapid for the case of periodic transverse 
boundary conditions than free transverse boundary conditions. This is similar to what was 
found in [31] and is due to the fact that periodic transverse boundary conditions minimize 
finite-size artifacts. Quantitatively, for a moderate value of the width, such as L y = 6, and 
q ~ 8, W is already within a few parts per thousand of its infinite 2D lattice value. 

TABLE III. Values of W for infinite-length, finite-width strips of the sqd lattice with periodic transverse 
boundary conditions, as functions of q, from evaluation of our exact solutions. For each pair of values of L y 
and q, the upper entry is the value of W from the exact solution and the lower entry is the ratio Rwp- 
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Ly [ q - 


6 


7 


8 


9 


10 


3 


1.817 


2.8845 


3.915 


4.932 


5.944 




1 


1 


1 


1 


1 


4 


2.629 


3.5005 


4.412 


5.348 


6.300 




1.024 


1.014 


1.009 


1.006 


1.004 



TABLE IV. Values of W for infinite-length, width L y strips of the sqd lattice with (PBC y , BC X ), as 
functions of q, compared with the corresponding values for the infinite 2D sqd lattice. For each value of L y 
and q, the upper entry is the value of W and the lower entry is the ratio of this value divided by W(sqd, q) 
for the 2D sqd lattice. For L y — 3, 4, the values of W are taken from the exact solutions; for 5 < L y < 10 
the values are from (6.4) with (6.5). 



L y [ q 


6 


7 


8 


9 


10 


3 


1.82 


2.88 


3.91 


4.93 


5.94 




0.743 


0.857 


0.908 


0.936 


0.953 


4 


2.63 


3.50 


4.41 


5.35 


6.30 




1.07 


1.04 


1.02 


1.015 


1.01 


5 


2.32 


3.29 


4.26 


5.23 


6.21 




0.949 


0.978 


0.989 


0.994 


0.996 


6 


2.43 


3.35 


4.29 


5.25 


6.22 




0.994 


0.994 


0.996 


0.997 


0.998 


7 


2.39 


3.33 


4.28 


5.25 


6.22 




0.976 


0.989 


0.994 


0.996 


0.998 


8 


2.41 


3.33 


4.29 


5.25 


6.22 




0.984 


0.991 


0.995 


0.997 


0.998 


9 


2.40 


3.33 


4.29 


5.25 


6.22 




0.980 


0.990 


0.994 


0.997 


0.998 


10 


2.40 


3.33 


4.29 


5.25 


6.22 




0.982 


0.990 


0.995 


0.997 


0.998 


oo 


2.45 


3.37 


4.31 


5.27 


6.24 




1 


1 


1 


1 


1 



One can also compare the results for W with those on the square and triangular lattices; 
this comparison shows the effect of the addition of a single diagonal bond to each square 
to go from the square to the triangular lattice, and then the addition of a second opposite 
diagonal bond to each square to go from the triangular to sqd lattice. This was done before 
in 
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VII. ZERO-TEMPERATURE CRITICAL POINTS 



TABLE V. Properties of P, W, and B for strip graphs G s of the sqd lattice. The properties apply for 
a given strip of type G s of size L y x L x ; some apply for arbitrary L x , such as N\, while others apply for 
the infinite-length limit, such as the properties of the locus B. For the boundary conditions in the y and x 
directions (BC y , BC X ), F, P, and T denote free, periodic, and orientation-reversed (twisted) periodic, and 
the notation (T)P means that the results apply for cither periodic or orientation-reversed periodic. The 
column denoted eqs. describes the numbers and degrees of the algebraic equations giving the XG s .j] f° r 
example, {6(1), 2(2), 2(3)} means that there are six linear equations, two quadratic equations and two cubic 
equations. The column denoted BCR lists the points at which B crosses the real q axis; the largest of these 
is q c for the given family G s . The notation "none" in this column indicates that B does not cross the real q 
axis. Column labelled "SN" refers to whether B has support for negative Re(q), indicated as yes (y) or no 
(n). 



Ly 


BCy 


BC X 


N x 


eqs. 


BCR 


SN 


1 


F 


F 


1 


{1(1)1 


none 


n 


2 


F 


F 


1 


{1(1)1 


none 


n 


3 


F 


F 


2 


{1(2)} 


none 


n 


4 


F 


F 


4 


{1(4)} 


none 


n 


1 


F 


(T)P 


4 


{2(1),1(2)} 


0, 2, 3 


n 


2 


F 


(T)P 


Q 

o 


{3(1)} 


0, 2, 4 


n 


3 


F 


(T)P 


16 


{6(1),2(2),2(3)} 


0, 2, 4, 4.25 


n 


3 


P 


F 


1 


{1(1)} 


none 


n 


4 


P 


F 


3 


{1(3)} 


none 


n 


3 


P 


(T)P 


4 


{4(1)} 


0, 2, 4, 6 


n 
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We summarize some of the results obtained in this work in Table [V]. Note that the 
L y = 1 strips of the sqd lattice are equivalent to L y = 2 strips of the triangular lattice, as 
discussed above. At the value of q where the nonanalytic locus B crosses the positive real 
axis, the g-state Potts antiferromagnet has a zero-temperature critical point [ 41] , |4| , [|5 1 . From 
the exact solutions for the chromatic polynomials and W functions, we thus conclude that 
the q = 2 Potts (Ising) antiferromagnet has a T = critical point on the L x — > oo limits 
of the cyclic/Mobius strips with L y = 2 and L y = 3 as well as the L y = 3 strip with torus 
(equivalently Klein bottle) boundary conditions. Note that this involves frustration owing to 
the non-bipartite nature of these graphs. As is generic for a T = critical point, the critical 
singularities are essential, rather than algebraic, singularities, as we have verified from an 
explicit transfer matrix calculation. The use of strip graphs with periodic or twisted periodic 
longitudinal boundary conditions is useful since the crossings of B on the real q axis signal 
the presence of T = critical points for the Potts antiferromagnet at these values of q. Just 
as was the case with the square and triangular strips, for which the Ising antiferromagnet 
also has a T = critical point |35| , ^ll , ^4| , ^5| , if one uses free longitudinal boundary conditions, 
B does not, in general, cross the positive real axis at q = 2. This difference is associated 
with the noncommutativity in the definition of W, as discussed before []12| , ^T1 , ^ , ^5| . We 
also find, for both the cyclic/Mobius strips and the torus/Klein bottle strips for which we 
have calculated exact solutions for W, in the nondegenerate cases, L y > 2, that the g-state 
Potts antiferromagnet has a T = critical point at q = 4. Since the L x — > oo limit of 
the cyclic/Mobius strips can be carried out with increasing even values of L x , for which the 
chromatic number is 4, this zero-temperature critical point is unfrustrated for these strips. 
In contrast, for the L x — > oo limit of the torus/Klein bottle strip, since min(x) = 6, it is 
frustrated. Finally, there is formally a similar critical point at q — 0. 



VIII. CONCLUSIONS 



In this paper we have presented exact calculations of the zero-temperature partition 
function (chromatic polynomial) and W(q), the exponent of the ground-state entropy, for 
the g-state Potts antiferromagnet with next-nearest-neighbor spin-spin couplings on strips 
of the square lattice strips (equivalently, the nearest-neighbor Potts model on strips of the 
sqd lattice) with width L y = 3 L y = 4 vertices and arbitrarily great length L x vertices. Both 
free and periodic boundary conditions are considered. In the L x — > oo limit, the resultant W 
function was calculated. By comparing the values of the exact W functions thus obtained 
for strips with various widths and boundary conditions versus numerical measurements of W 
for the full 2D sqd lattice, we evaluated the effects of the next-neighbor spin-spin couplings. 
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We showed that the q = 2 (Ising) and q = 4 Potts antiferromagnets have zero-temperature 
critical points on the L x — > oo limits of the strips that we studied. With the generalization 
of q from Z + to C, we also determined the analytic structure of W(q) in the q plane. 
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IX. APPENDIX 
A. A Reduction Theorem 

Consider a strip of a given width and length, made up of squares such that the four 
vertices of each square are connected to form a complete graph K r (so that there are r — 4 
vertices outside of the strip) with some set of boundary conditions (BC). We shall denote this 
strip temporarily as G sq ^K r ,BC- We first show that the chromatic polynomial for this graph 
can easily be expressed in terms of the chromatic polynomial for the corresponding graph 
with the vertices of each square forming a K±, i.e. such that the two diagonally opposite 
vertices of each square are connected by an edge. For this purpose we use the intersection 
theorem from graph theory; this states that a graph G can be expressed as the union of two 
subgraphs G = G\ U G2 such that the intersection of these subgraphs is a complete graph, 
G\ nG 2 = Kj for some j, then P(G, q) = P(Gi, q)P(G 2 , q)/P(Kj, q). Applying this theorem 
to each square of the above strip, we have 

P(G sq , Kr , BC , q) = ( J P(G sg , K4 , BC , q) = " j)\ P(G sg ,K 4 ,BC, q) (9.1) 

where N4 denotes the number of squares on the strip and we use the standard notation from 
combinatorics for the "falling factorial" , 
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s-1 

ff W = JJ(g - j) . (9.2) 

j=0 

Our results in the text thus also apply to lattices comprised of squares such that the four 
vertices of each square form a K r with r > 4. 

B. (K r , K s ) Strips 

We discuss here a strip of t K r subgraphs, connected such that successive K r subgraphs 
intersect on a K s subgraph, with 1 < s < r — 1, with some longitudinal boundary conditions 
imposed. A member of this general family may be labelled as (K r , K s ,t, BC X ). Two of 
the simplest longitudinal boundary conditions to impose are free and cyclic, which we shall 
denote as FBC X and PBC X . The general (K r , K s , t = m, FBC X ) graph has n = (r — s)m + s 
vertices. 

One simple set of families is the cyclic strip (K r , Ki, t = m, PBC X ). These may think of 
these graphs as being formed by starting with a circuit graph, C m and gluing to each edge 
one edge of a complete graph K r . For these we find 

P{{K r ,K x ,m,PBC x ),q) = (^jP(C m ,q) = [f[(q - j)]™ P(C m ,q) (9.3) 

J— ^ 

where P(C m , q) is the chromatic polynomial for the circuit graph with m vertices, given in 
the introduction. 

Another interesting family is the open strip (K r , K s ,m, FBC X ). We calculate 

P((K r ,K s ,m,FBC x ),q) = q^(^j = g « [JJ( g - j)\ . (9.4) 

/ j =s 



Hence, in the m — ► oo limit, 

q 



w 



(r)N — r-1 



[n^-^)] r_s ( 9 - 5 ) 



J=S 



and B = 0. 

Other cases are more complicated. For a strip of -KVs, i.e., triangles, in addition to 
free and periodic (cyclic) boundary conditions, one also has the possibility of twisted peri- 
odic boundary conditions, which form a Mobius strip. The chromatic polynomial for the 
(K 3 , K 2 , t, PBC X ) strip for even t was given in (see also f3~8fl); in this case, the strip can 
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be regarded as being formed from a strip of m squares with a bond added to each square 
connecting the lower left to upper right vertex of the square, so that t = 2m. The chro- 
matic polynomial for the corresponding Mobius strip, (K 3 , K 2 ,t,TPBC x ) with even t was 
also given in [j35] . For the t — > 00 limit, the W function and associated nonanalytic locus B 
were given in |35|]. The chromatic polynomials for the corresponding strip with odd t, viz., 
(K3, K 2 ,t, PBC X ) and (K3, K 2 ,t,TPBC x ), were given in f45fl ; the t — > 00 function W and 
locus B are the same for both the cyclic and Mobius strips as t increases through even or 
odd values. 

In the text we have considered strips in which each square forms a K4 subgraph. There are 
actually different classes of strips with t repeated subgraphs, in the notation introduced 
above. An illustration of these is given in Fig. |8|. 




FIG. 8. Illustrative strip graph comprised of K4 subgraphs intersecting on common edges (-KVs), of a 
type different than that shown in Fig. 1(a). 



As an example, we calculate 

P((# 4 , K 2 , t, BC X ), q) = (q- 3)* P((K 3 , K 2 , t, BC X ), q) 



(9.6) 



where BC X = PBC X or TPBC X , respectively and the (K 3 , K 2 ,t, PBC X ) and 
(K 3 , K 2 ,t,TPBC x ) strips are the cyclic and Mobius triangular-lattice strips mentioned 
above. Thus, explicitly, for even t = 2m, 

P({K*, K 2 , t = 2m, PBC X ), q) = {q - 3) 2m [(q 2 - 3g + 1) + [(q - 2) 2 } m + {q - 1) [(\ t2 , 3 ) m + (A t2j4 

(9.7) 



P((K 4} K 2 , t = 2m, TPBC X ) } q) = {q- 3) 2m -1 + [(q - 2) 2 } m - (q - l)(q - 3) 



where X t2 j, j = 3,4, were defined in eqs. (|1.7|) above. 
For the case of odd t = 2m + 1, we have 



A*2,3 — At2,4 

(9.8) 
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P((K 4 ,K 2 ,t = 2m + l,PBC x ),q) = (q - 3) 



2m+l 



V-3g+l) + (g-2)[(g-2)r+ 



:(?-i)(9 -3) 



(A t2>3 ) m + (A t2 , 4 ) m ) 



+ 



(Ai2,3) m — (^2,4) 
A<2,3 — At2,4 



(9.9) 



P((K 4 , K 2 , t = 2m + l, TPBC X ), q) = (q - 3) 



2m+l 



l + (g-2)[(g-2) 2 r+ 



((V 3 ) m + (Ai 2 , 4 ) m ) + (9 - 4g) 



(At2,3) m — (A<2,4) r 
At2,3 — Aj2,4 



(9.10) 
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